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Recent approaches for quantum gravity are conjectured to give predictions for a minimum measur- 
able length, a maximum observable momentum and an essential generalization for the Heisenberg 
uncertainty principle (GUP). The latter is based on a momentum-dependent modification in the 
standard dispersion relation and leads to Lorentz invariance violation (LIV). The main features of 
the controversial OPERA measurements on the faster-than-light muon neutrino anomaly are used 
O | to calculate the time of flight delays At and the relative change Av in the speed of neutrino in 

dependence on the redshift z. The results are compared with the OPERA measurements. We find 
that the measurements are too large to be interpreted as LIV. Depending on the rest mass, the 
propagation of high-energy muon neutrino can be superluminal. The comparison with the ultra 
high energy cosmic rays seems to reveals an essential ingredient of the approach combining string 



a 

CD 

' theory, loop quantum gravity, black hole physics and doubly spacial relativity and the one assuming 

, a pertubative departure from exact Lorentz invariance. 
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ON I The combination of quantum mechanics and special relativity reveals serious difficulties in describing the 

one-particle theories. As combining the Hciscnbergs uncertainty principle and finiteness of the speed of light c 
leads to creation and annihilation processes, when studying the Compton wavelength of the particle of interest 
Q. Another consequence of the space-time foamy structure at small scales is the Lorentz invariance violation 
(LIV). Different approaches for the quantum gravity Q, the yet-to-be-built quantum theory of gravity, have 
been proposed [H, y]. They provide a set of predictions for a minimum measurable length, a maximum ob- 
servable momentum and an essential modification of the Heisenberg uncertainty principle (GUP). The latter is 
proportional to a quadratic momenta according to string theory, loop quantum gravity and black hole physics. 
Based on doubly spacial relativity, a proportionality to first order moments (linear) has been suggested. Ali's 
model [H, which is utilized in this work, combines both approaches. 

The suggestion that Lorentz invariance (LI) may represent an approximate symmetry of nature dates back 
to about four decades @ . This represents the roots of LIV. A self-consistent framework for analyzing possible 
violation of LI was suggested by Coleman and Glashow [3, 0. On a sample of gamma ray bursts (GRB), the 
energy dependent time offsets are investigated in different energy bands assuming standard cosmological model 
. A weak indication for redshift dependence of time delays suggestive of LIV has been found. A comprehensive 
review on the main theoretical motivations and observational constraints on Planck scale suppressed Lorentz 
invariance violation is given in [l(| and the references therein. The Planck scale itself turns to be accessible in 
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Following the proposal to utilize astrophysical objects to look for energy-dependent time of arrival delays 
we present in the present work an estimation for the time of flight delays and the relative change in the 
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velocity of muon neutrino based on GUP-approach which will be introduced in section [TTJ Section Hill is devoted 
to the calculations which are confronted with the recent OPERA measurements and ultra high energy cosmic 
rays (UHECR) observations. The conclusions are addressed in section HVl 



II. THE APPROACH 



For a particle with mass M having a distant origin and an energy scale comparable to the Planck's one, 
the momentum would be a subject of a tiny modification [1, According to GUP-approach, the comoving 
momenta can be given as 



p v (l - ap + 2a 2 po) , 



(1) 



where po is the momentum at low energy. The parameter a — a>o/(cM p i) = a^l v \jh [4, 5], where c, ao, M p i (l p i) 
are speed of light as introduced by Lorentz and implemented in special relatively, dimensionless parameter of 
order one, and Planck mass (length), respectively. Then, the dispersion relation in comoving frame reads 



E, 



pi c 2 (l-2ap ) + M 2 c 4 



(2) 



When taking into consideration a linear dependence of p on a and ignoring the higher orders of a, then the 
Hamiltonian is 



H = (plc 2 -2aplc 2 + M 2 u c A Y l . 
The derivative of H with respect to the momentum results in a comoving time-dependent velocity 



v(t) 
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(4) 



The comoving momentum can be related to the physical one, p v = p Uo (to)/a(t), where a is the scale factor, 
which in turn can be related to the redshift z 



a(z) 



1 



1 + 2 

In the relativistic limit, p » M, the fourth and fifth terms in Eq. ((4]) simply cancel each other. Then 

M, 2 c 2 Mtc A 



v{z) = c(l + z) 



l-2a{l + z)p v 



(5) 



(0) 



2(1 + z) 2 p 2 „ 2{l + zfpl o _ 
where po is treated as a comoving momentum. It is straightforward to deduce the relative change in the velocity 



Av(z) 



-2(l + z) 2 p„ 
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(7) 



The results are given in the left panel of Fig. [TJ For a masslcss muon neutrino, the sign of Av (z)/c remains 
negative with increasing z. When taking into account the muon neutrino mass, the sign turns to positive. In 
this case, its value nearly vanishes at large z. Accordingly, the resulting sign of the summation of first two terms 
of Eq. (J7|) is determined by the second term, i. e., positive, at small z. Then, the sign is flipped to negative at 
z ~ 0.2, i. e., the first term becomes dominant. 

The comoving redshift-dependent distance travelled by the particle of interest is defined as 



r{z) 



v(z) 



-dz. 



I (l + z)H(z) 

where H(z) is the Hubble parameter in dependence on z. From Eqs. ^ and (|5]l. the time of flight reads 



I -2a(l + z)p p 



M 2 c 2 



M*c 4 



2(1 + z) 2 Pl 
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dz 



(8) 
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which inheres the well-known time of flight of a prompt low-energetic photon (first term), which apparently is 
invariant in Lorentz symmetry. It is apparent that Eq. §§§ contains a time of flight delay 



2a (1 + z)p Vo 



4„4 



4(l + 2) 3 p3 



2 „2 



2(1 + ^)X 



dz 



(10) 



due to various LIV effects (first and second terms) and the particle mass (third term). The second term seems 
to reflect a mixture of LIV and the effects of the rest mass on delay time. 

In order to determine Eq. (TIT)]) , it is obvious that the redshift-dependence of the Hubble parameter H(z) 
plays an essential role. Furthermore, H(z) depends on a time-dependent redshift, dz/dt. 



da{z) 

H(z) = 



dz 
dt 



1 dz 



a(z) 



1 + z dt' 



(11) 



It is obvious that this expression can be deduced from Eq. ([5]). In general, the expansion rate of the Universe 
varies with the cosmological time fl3M20l ] . It depends on the background matter/radiation and its dynamics 
flfjj . The cosm olog ical constant reflecting among others the dark matter content seems to affect the temporal 



evolution of H j 1 81 ] . Fortunately, the redshift z itself can be measured with a high accuracy through measuring 
the spectroscopic redshifts of galaxies having certain uncertainties {a z < 0.001). Based on this, a differential 
measurement of time, dt, at a given redshift interval automatically provides a direct and clean measurement of 
H(z) [2114231 ] . These measurements can be used to derive constraints on essential cosmological parameters 24 1. 



In this letter, we plan to implement the measurements of the expansion rate and their constrains in evaluating 
the integrals given in Eq. (fTO")) . 

Out of a large sample of early- type galaxies (about 11000) extracted from several spectroscopic surveys 
spanning over ~ 8 x 10 9 years of cosmic lookback time, i. e., 0.15 < z < 1.42 22[, most massive, red elliptical 
galaxies, passively evolving and without signature of ongoing star formation are picked up and used as standard 
cosmic chronometers [24[. Then, the differential age evolution turns to be accessible. It gives a cosmic time 
and can directly probe H(z). A list of new measurements of H(z) with 5 — 12% uncertainty is introduced in 
[22I ], The uncertainty seems to be comparable with the one achieved for H(z ~ 0.2). Figure Q] illustrates these 
measurements in BC03 [2o| model and in combination with the cosmic microwave background (CMB) data and 
constraining possible deviations from the standard (minimal) flat ACDM model [23[. The right panel shows a 
data set taken from MS model (26|. It gives a clear evidence that the results are model-dependent. 
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Fig. 1: Left panel: the Hubble parameter H calculated in BC03 model (open triangle) and in combination with CMB 
data and constraining possible deviations from the standard (minimal) flat ACDM model (solid circles) is given in 
dependence on the redshift z. The results from MS model are drawn in the right panel. The curves represent the fitting 
parameters (see text for details). 

The measurements are fitted as follows. For the results obtained from BC03 model [25[ and using a com- 
bination with CMB data and constraining possible deviations from the standard (minimal) flat ACDM model 



H(z) =01+1/1* + h z 2 . 



(12) 
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where Pi = 72.68 ± 3.03, 71 = 19.14 ± 5.4 and Si = 29.71 ± 6.44. This expression is given by the solid curve 
(left panel of Fig. [T]). For the MS model (2(| measurements, we suggest two expressions 



H(z) 



P2+l2Z + 5 2 Z 2 +€ 2 Z 3 , 

03 + 73tanh(5 3 2:), 



(13) 
(14) 



where p 2 = 66.78 ± 8.19, 72 = 113.27 ± 7.5, S 2 = -140.72 ± 12.6, e 2 = 60.61 ± 5.48, f3 3 = 71.94 ± 4.35, 
73 = 33.51 ± 7.94 and £3 = 1.6 ± 0.1. The results of Eq. (fT3)) are given by dashed curve in the right panel 
of Fig. [1] Equation (fT4")l is drawn by dotted curve, where the largest point is excluded while other points 
are constructing the ensemble used in the fitting. It is obvious that the implementation of Eq. (|13|) . which is 
obviously a rational function, in Eq. (|10[) results is a non-analytic integral. Implementing of Eq. (|14[) in Eq. 
(II OP makes the second and third integrals non-solvable, while the first term cab be solved as given in Appendix 
lAl where the results are also given. 

It is apparent that Eq. (jT2"j) simplifies the integrals in Eq. (jTU)) . Accordingly, there are two types of LIV 
contributions to the time of flight delay. The first type is originated in LIV, i.e., finite a. It appears in two 
terms. 
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where A = (4/3i<5i — 7 2 ) 1 / 2 . Furthermore, Eq. (p~2]) results in an exclusive estimation for the mass contribution 
to the time of flight delay, 
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The results arc discussed in section PlI Al 



III. CONFRONTING WITH MEASUREMENTS AND OBSERVATIONS 



A. OPERA Measurements for At and Av 



The results of Eqs. (fT5j) and (|T6)) are given in the left panel of Fig. [2] In these calculations, we use 
the same configurations of the controversial OPERA experiment [27| in measuring a faster-than-light muon 
neutrino anomaly, i.e., M v = leV and the beam energy E v = 17GeV. Assuming that this neutrino beam 
is stemming from a distant origin and was witnessing such a huge redshift z, the time of flight delay At is 
calculated with varying z. The first two terms of Eq. (jTUJ) are calculated and drawn, separately. We find that 
the first term 0(a) is one or two orders of magnitude higher than the second one 0(a,M). Their summation 
can be approximated as At ~ 10~ 12 sec. For a neutrino beam covering the distance between the source of 
CERN neutrino beam directed to Gran Sasso (CNGS) and the OPERA detector at the underground Gran 
Sasso Laboratory (LNGS), ~ 730km, a time of flight delay of ~ 1045.1 ± 11.3 nano seconds is believed to be 
registered. Despite of its statistical interpretation, when comparing this with the LIV time of flight delay, we 
find that OPERA measurement is too huge (six orders of magnitude) to be understood as LIV. 

Furthermore, the OPERA experiment suggests an increase in the speed of light by about 7.5 km/sec (~ 25 
part per million of c) [27} • The left panel of Fig [2] illustrates the redshift evolution of the possible change in the 
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Fig. 2: Using the main features of OPERA experiment, the relative change in the velocity of muon neutrino is given as 
a function of redshift z in the left panel. The right panel shows the time of flight delay. 



velocity of muon neutrino having mass M = 1 eV and energy 17 GeV according to LIV. The first two terms of 
Eq. ([7]) are compared with each other, while z raises. We notice that the first term (massless muon neutrino) 
has a negative sign. Its absolute value increases almost linearly with increasing z. It is obvious that the second 
term remains in the positive site of the ordinate. Its value decays almost exponentially, while abscissa raises. 
The upper and lower values are ranging from ~ 10~ n to 0. In order words, the second term expresses a positive 
change in the neutrino velocity. Quantitatively, the change should be less than ~ 10~ n c. 

With increasing z, the sum of these two terms results in a flipping in the sign of Av(z)/c. At small z, the 
second term seems to be dominant. At z ~ 0.2, the positive sign is switched into negative. At larger z-valucs, 
the first term becomes dominant. So far, we conclude that the average change in the velocity of muon neutrino 
is Av ~ —2 x 10~ n c, which is almost six orders of magnitude smaller than the value measured in OPERA [27j |, 
That the sign of Av(z) = c v — c is flipped would mean that 

{C(+ivc), then c v = c + O 
(18) 
O(-ivc), then c„ = c-O 

where c v is the velocity of neutrino. From Eq. (|18p . the second case apparently follows the Lorentz invariance 
symmetry. The first case suggests that the speed of light would not be constant in all inertia frames. Further- 
more, it would not be the maximum of travelling matter and information in the universe. The value of O is 
about ~ 10 -11 c, i. e., ~ 3 mm/sec. Regardless its small value, it seems to describe a superluminal propagation 
of high-energy muon neutrino at z < 0.2. 

It is believed that OPERA is giving results comparable with the MINOS [28| earlier ones. The values of the 
velocity change, O, have been assumed as ~ 10 -5 . On the other hand, these measurements are not compatible 
with the observations of ~ 10 MeV-neutrino from supernova SN1987a [2!|, where the value of O is measured 
as ~ 10 -9 . Therefore, the faster-than-light anomaly seems to be energy-dependent. It drops rapidly, when 
reducing energy from GeV- to MeV-scale Q. Nevertheless, the velocity anomaly is conjectured to reflect the 
propagation of all decay channels of neutrino and new physics such as LIV. With this discussion, we want to 
emphasize that the recent OPERA measurements might far away to construct a solid reference. 



B. Ultra High Energy Cosmic Rays 

Other consequences of LIV can be tested in ultra high energy cosmic rays (UHECR). Following from Eqs. 
© and ([TU]) . the time of flight is delayed by the factor At, then the pair production is kinematically allowed, 
when [| 

/ 2 \ 1/2 

E 7 >m e ( — \ , (19) 
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where the subscript e stands for electron or positron and V = (v e — c)/ c. According to Eq. (fT8|) and assuming 
that the electron has the same characteristics as the muon neutrino of OPERA experiment. The electron velocity 
can be estimated as (1+1?) c, where T> can be ne gati ve or positive. That photons with energies E 1 > 50 x 10 12 eV 
[30| and electrons with energies £ e ~lx 10 12 eV [3lj have been observed would set upper limits to V, < 2 x 10 -16 
and T> ~ 1.3 x 10~ 13 , respectively. It is apparent that these values are energy-dependent. Both are smaller 
than Av ~ 10~ n c estimated by GUP. Such a discrepancy would be interpreted as follows. In our calculations 
we apply a GUP-approach depending on linear modification on momentum 0, [El- As discussed above, this 
approach combine string theory, loop quantum gravity, black hole physics and doubly spacial relativity. Stecker 
and Glashow use 7-ray observations of Mrk 501 in constraining the Lorentz invariance breaking parameter || 
based on potential departure from exact Lorentz invariance introduced in a perturbative framework Q . 

In light of recent theoretical work on quantum gravity, especially within string theory, the sensitivity factor 
of gamma ray bursts (GRB) reads 77 = \At*\/5t [12j, where St reflects the time structure of the signal and At* 
is a special time of flight delay, when E qg ss E p i, where E qg (E p i) is an effective quantum gravity energy scale 
(Planck energy scale). It is obvious that At* would be equivalent to At of this work. Depending on distant origin, 
GRB emission can reach the Earth with different time structures St. The conventional gravitational lensing 
is achromatic. Therefore, the energy- dependent time delay would not be dependent on the actual emission 
mechanism GRB. A lensed GRB observed couple decades ago [HJ cab be used to estimate the sensitivity factor 
as 77 « 10~ 6 revealing that St rj 10 _7 scc. When pulsars, supernovae and other astrophysical phenomena, but 
not GRB, are considered (HI, 77 « 10 -10 and therefore St « 10 _3 scc. For neutrinos from type-II supernovae 
like SN1987a, 77 w 10~ 4 and the time structure can be estimated as « 10 -9 sec. In principle, the upper bound 
on a parameter by comparing the calculations with the experiments We would do this with OPERA and/ 
or MINOS. GRB would allow us to set an upper value for the GUP-charactering parameter a (Hf. 



IV. CONCLUSION AND OUTLOOK 



In this paper we performed calculations for the time of flight delays and the relative change in the velocity 
of muon neutrino with mass leV and energy 17GeV. In doing this, we utilized the GUP-approach, which is 
based on a momentum-dependent modification in the standard dispersion relation. For a particle having a 
distant origin and energy comparable with the Planck energy scale, the comoving momentum is given as a 
series of linear modifications on momentum. We calculated the relative change in the speed of massive muon 
neutrino and its time of flight delays with varying redshift. The latter depends on the temporal evolution 
of the Hubble parameter, which can be estimated from a large sample of early-type galaxies extracted from 
several spectroscopic surveys spanning over ~ 8 x 10 9 years of cosmic lookback time, most massive, red elliptical 
galaxies, passively evolving and without signature of ongoing star formation are picked up and used as standard 
cosmic chronometers giving a cosmic time directly probe for H(z). The measurements according to BC03 
model and in combination with CMB data and constraining possible deviations from the standard (minimal) 
flat ACDM model are used to estimate the z-dependence of Hubble parameter. The measurements based on 
MS model are used to show that the results are mo del- dependent. 

We compared the results with the OPERA experiment. Thus, we conclude that the OPERA measurements 
for At and Av arc too large to be interpreted as LIV. Depending on the rest masses, the propagation of high- 
energy muon neutrino can be superluminal. The comparison with UHECR reveals the potential discrepancy 
between an approach combining string theory, loop quantum gravity, black hole physics and doubly spacial 
relativity and a pertubative departure from exact Lorentz invariance. 

In light of this study, we believe that GRB would be able to set an upper value for the GUP-charactering 
parameter a. Furthermore, the velocity anomaly is conjectured to reflect the propagation of all decay channels 
of neutrino and new physics such as LIV. 
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Fig. 3: Left panel: the time of flight delay of muon neutrino (mass 1 eV and energy 17 GeV) in dependence on z, using 
(|13|l in Eq. (|10|) . The first term is given by dashed curve, while dotted curve represents the second term. The right 
panel draws Eq. (|A1|) . Both dashed curves seem to represent comparable results. 



Appendix A: Time of Flight Delay According to MS Model 



It is apparent that integrating the rational expression (|13[) into Eq. (|10p gives a numerical solution. In left 
panel of Fig. [3l the first (dashed curve) and second (dotted curve) terms of Eq. (fTO)) . where H (z) is taken from 
(fT3"l) . are given in dependence on z. Their summation is given by the solid curve. The time of flight delay, At 
can be averaged as ~ 10 _13 sec. This value is much smaller than the one measured in OPERA experiment, so 
that the latter would not interpreted by LIV. 

When implementing Eq. (|14[) into Eq. (fT0|) . the integrals in the second and third terms can not be solved, 
analytically. The first term can be solved as follows. 



At(z) = 2a Pvo 



.„ . dz 
H(z) 

ap VQ B sech (S3 z) 



PlSle c z(2 + z) + 



e c fa (Pi - 7!) <5f [ft + 73tanh (S3 z)] 
(3 3l3 e c (iirln[l + e 2S3Z ] -2Cln(l- 
2 5 3 [^z + zC + z\n(l-e^ c+5 >^ 



„S 3 z-2C\ 



i 7r In [cosh (S; 




-e C hiStz* (Al) 
2Cln [isiah(0 + S 3 z)} 



]nB 



})]} 



where B = /?3Cosh (S3 z) + 73sinh (S3 z) and C = atanh (fa/ ^3). The results arc drawn in the right panel of Fig. 
[31 In these calculations, only the real component of the second line of Eq. (|A1[) is taken into consideration. 
The values of At can be approximated to ~ 10~ 13 sec, which about seven orders of magnitude smaller than the 
time of flight delay measured in OPERA experiment. With the dashed curve (first term) in left panel, this term 
gives comparable results, qualitatively and almost quantitatively. 
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